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1.  INTRODUCTION 

Julia  sets  are  amongst  the  most  frequently  pictured  fractals,  combining  both  aesthetic  and  actual  beauty.  The 
modem  day  interest  in  Julia  sets  and  related  mathematics  began  in  the  year  1920,  which  was  initiated  by  French 
mathematician  Gaston  Julia.  In  1918  he  wrote  a paper  titled  "M'emoire  sur  l’ iteration  des  functions  rationelles"  (A  Note 
on  the  Iteration  of  Rational  Functions)  [11]  where  he  first  introduced  the  modem  idea  of  a Julia  set.  In  this  paper  Julia  gave 
a precise  description  of  the  set  of  those  points  of  the  complex  plane  whose  orbits  under  the  iteration  of  a rational  function 
stayed  bounded.  Interest  in  the  subject  flourished  over  the  next  10  years  and  many  other  well-known  mathematicians  like 
Harald  Cramer  began  to  study  Julia  set.  Despite  the  lack  of  computing  machine  available  at  that  time  he  was  able  to 
become  the  first  man  to  approximate  an  image  of  a Julia  set.  Due  to  the  lack  of  computing  machine  the  progress  of  research 
in  this  line  slowed  down  and  after  some  years  Julia's  work  was  forgotten  by  the  mathematical  community.  It  was  the 
France  mathematician  Benoit  Mandelbrot  who  brought  back  Julia's  work  around  1977.  With  the  aid  of  computer  graphics 
he  showed  that  Julia's  work  is  a source  of  some  of  the  most  beautiful  fractals  known  today. 

A very  brief  historical  note  will  be  useful  to  explain  the  role  of  iteration  of  rational  maps  on  the  complex  plane 
which  is  the  central  theme  of  this  paper.  The  study  of  iteration  of  maps  begins  with  the  Newton's  method.  In  1870's  Cayley 
and  Schroder  independently  studied  Newton's  method  over  the  complex  plane  [7],  Julia's  interest  on  iteration  of  maps 
apparently  was  motivated  by  the  1879  paper  "The Newton-Fourier  Imaginary  Problem"  by  Sir  Arthur  Cayley  [14],  Cayley 

was  studying  the  equation  f (z)  = z3  + C = 0 , using  the  Newton-  Fourier  iterative  method  to  find  the  roots.  In  this 
iterative  method  one  can  approach  to  a root  as  a limiting  value  for  a certain  starting  value  of  z . Since  there  are  three  roots, 
Cayley  tried  to  find  out  a method  by  which  it  can  be  ascertained  to  which  of  these  three  roots  the  iteration  will  converge  for 
a given  starting  value  of  z . But  he  could  not  solve  this  problem.  Over  the  following  decades  other  mathematicians 
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continued  to  build  on  this  work,  but  their  work  failed  to  create  a general  and  global  theory  that  could  describe  the  behavior 
of  all  rational  functions  on  the  entire  complex  plane. 

In  1915  the  French  Academy  of  Sciences  announced  that  it  would  award  the  1918  Grand  Prix  des  Sciences 
Mathematique  for  such  a general  and  global  theory.  The  two  mathematicians  who  succeeded  in  creating  the  modem  theory 
of  complex  dynamics  were  Gaston  Julia  and  Pierre  Fatou.  Julia  won  the  contest  and  Fatou,  despite  withdrawing  his  work 
from  the  contest,  received  a second  prize.  Both  saw  that  iteration  of  any  function  partitioned  the  extended  complex  plane 
into  a region  where  the  iterates  are  equicontinuous  and  a region  where  they  are  not.  Today,  for  any  non-constant  rational 
function  f , the  maximal  open  subsets  of  the  extended  complex  plane  on  which  the  family  of  iterates  of  the  function  f , 

i.e.,  \f" : n > l}  is  equicontinuous  is  called  Fatou  set  of  f , denoted  by  F or  F (/)  , while  the  complement  of  F in 
the  extended  complex  plane  is  called  the  Julia  set  of  f , denoted  by  J or  J (f) . This  concept  is  the  primary  basis  of  this 
paper. 

With  the  arrival  of  modem  sophisticated  computers  in  1980s  the  above  mentioned  intricate  problem  of  course 
could  be  handled  and  the  research  in  this  area  got  momentum.  Let  us  see  little  bit  in  detail  the  mathematical  theory  behind 
the  above  mentioned  problem. 

f(z) 

Consider  F ( z ) = z - - — \ . 

f\z ) 

Suppose  z0  be  a root  of /(  z ) , i.e.  / ( z0  ) = 0 . Then  we  have  F ( z0  ) = z0  , which  shows  that  z0  is  a fixed 
point  ofF(z  ). 

Thus,  any  complex  number  z is  a root  of  / ( z ) = 0 if  and  only  if  it  is  a fixed  point  ofF  ( z ) . 

/ \ z3  - 1 2z3  - 1 

Now,  if  we  consider  F(z  ) = z — = — . 

3z2  3 z2 

2 ni  Ani 

Then  clearly  the  roots  of  /(z)  are  the  cube  roots  of  unity,  viz.,zj  =1,  z2  = e 3 and  z3  = e 3 . These  are 
also  attracting  fixed  point  ofF  ( z ) . Therefore,  the  problem  of  interest  is  now  converted  to  find  out  the  basin  of  attraction 
of  these  fixed  points  of  F ( z ) . Using  a computer  to  graph  this  situation  in  the  complex  plane  we  get  the  image  as  shown  in 
the  figure  1.1.  Here,  the  basin  of  attractions  forz1;  z2  and  z3  are  colored  respectively  by  blue,  green  and  yellow.  Note 
that  the  basin  boundary  in  this  figure  is  nothing  but  the  Julia  set  for  the  function  F ( z ) . 
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Figure  1.2 


With  a similar  fashion  one  can  find  the  basin  of  attraction  for  the  roots  of  the  equation  z4  - 1 = 0 as  shown  in 
the  figure  1.2,  where  the  region  colored  in  blue,  yellow,  green  and  red  are  respective  basin  of  attraction  for  the  roots  of 
/ ( z ) or  fixed  points  of  f(z)  which  are  1,-1,  i and  - i . 

For  fundamental  work  relevant  to  study  of  iteration  of  rational  polynomials  one  may  refer  Beard  on  [1],  Devaney 
[5,  6]  and  Steinmetz  [18].  For  computing  Hausedorff  dimension  of  the  Julia  sets  one  can  see  [2],  Julia  set  is  a classical 
example  of  fractal,  which  has  lots  of  application  in  medical  science,  study  of  weather  system,  biology,  wireless 
communications  etc.  Interested  readers  can  go  through  [ 10,  12]  for  finding  such  applications.  Julia  set  which  is  a fractal  is  a 
great  tool  for  animating  the  various  aspect  of  world  around  us,  it  is  used  to  animate  feature  films  like  'Star  Track-II1  [12]. 
Moreover,  for  applications  of  Julia  sets  in  wavelet  analysis  and  in  the  study  the  large  scale  distribution  of  galaxies  in  the 
observed  universe  one  can  see  [13]  and  [16]  respectively. 

The  rest  of  the  paper  is  organized  as  follows.  In  section-2,  we  provide  a review  of  preliminary  concepts  and 
definitions  along  with  some  examples  to  clarify  those  concepts  and  definitions.  Section  -3  and  4 deals  with  some 
dynamical  and  structural  properties  related  to  Julia  set.  In  section  5 we  have  described  about  how  to  visualize  the  Julia  set. 
Finally,  we  have  given  our  conclusions  in  section  6. 

2.  SOME  PRELIMINARIES 

As  discussions  on  Julia  set  involve  the  extended  complex  plane  (that  is  infinity  is  included),  we  have  given  a brief 


www.iaset.us 


editor@iaset.us 


102 


Arun  Mahanta,  Hemanta  Kr.  Sarmah,  Ranu  Paul  & Gautam  Choudhury 


review  of  how  it  can  be  achieved. 

Stereographic  Projection:  Stereographic  projection  is  a method  to  create  a planar  map  of  a sphere.  In  the  figure 
given  below,  we  draw  a sphere  of  radius  one  whose  south  pole  touches  the  complex  plane  at  the  origin.  To  each  point  of 
the  sphere  there  is  a corresponding  point  in  the  complex  plane.  To  be  specific,  the  ray  from  the  North  Pole  through  the 
point  ( X,y,z ) on  the  sphere  will  intersect  the  complex  plane  at  a point  denoted  by  U + iv  enabling  us  to  define  a map  as 
2 

P(x,y,z)  = - (. x+iy ) 

2 — z 


Figure  2.1:  Stereographic  Projection 

2 2 

Hence,  u = x and  v = y . 

2—z  2-z 

The  inverse  of  this  transformation  is  given  by 

P~{  ( u + iv ) = — - — ^ (■ 4v,4w,2m 2 +2v2). 

u + v~  +4 

Which  gives  the  point  (x,  y,  z)  on  the  sphere  corresponding  to  the  point  It  +iv  of  the  complex  plane  C The  z- 
component  of  this  inverse  transformation  is 

z=4"!±^ 

u2  +v2  +4 


. .22 

When  U + IV  tends  to  infinity,  i.e.,  U + V goes  without  bounds  the  z-component  tends  to  2 since, 


2 (»2+v2) 


lim  — f r — 

a2+v2-X»  u~  4-  V +4 


= 2 But,  the  only  point  on  the  sphere  whose  z-component  is  2 is  the  North  Pole  i.e.,  the  point 


(0, 0, 2) . Hence,  the  North  Pole  can  be  regarded  as  the  point  at  infinity. 
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In  this  way  the  extended  complex  plane  can  be  considered  as  a sphere,  called  Riemannian  sphere.  We  denote  it  by 

C , i.e.,  C = C VJ  {co}. 

Let  <J  be  the  metric  on  the  Riemannian  sphere  which  simply  denote  the  length  of  the  chord  between  any  two 
points  on  the  sphere.  This  metric  is  called  spherical  metric  or  the  chordal  metric,  and  for  any  two  points  z , W G C , it  is 
defined  as 


a(z,w)  = 


2z-w 


V1+N2  -J 


1 + \w\ 


When  w — > go  , then  <7{z,w ) = 


1+  z' 


With  this  metric  C is  a compact  metric  space  [14], 

Consider  an  arbitrary  polynomial  f \ C — > C . Let  f " denote  the  nth  iterate  of  f , that  is  , f composed  with 
itself  n times.  For  each  point  z0  e C , we  are  interested  in  the  behavior  of  the  sequence 
Z0 , f{zQ  ),  f2  (zQ  ),•••,  f"  (z0  ),  ■ ■ ■ , and  in  particular,  what  happens  as  Tl  goes  to  infinity. 


Now  we  review  several  definitions  related  to  the  iterations  of  rational  functions  which  are  relevant  to  carry  our 


study. 


Definition  2.1:  A function  f is  called  analytic  in  C if  its  derivative  exists  at  each  point  in  C . 

Definition  2.2:  A point  z0  e C is  called  critical  point  for  the  function  f if  f\z0)  — 0. 

Definition  2.3:  A point  Z0eC  is  called  periodic  point  of  f if  f"  (zQ ) = z0  for  some  integers  > 1 . The 
smallest  n with  this  property  is  called  the  period  ofz0.  Thus,  the  periodic  points  of  z0  are  the  zeros  of  the  function 
F(z0,f)=fn(z0)-z0. 

A periodic  point  with  period  one  is  termed  as  fixed  point  of  f i.e.,  Z0  is  a fixed  point  of  / it  f (z0  ) = Z0  . 

Definition  2.4:  The  multiplier  (or  Eigen  value,  derivative)  A of  a rational  map  f iterated  U times,  at  the 
periodic  point  z0  is  defined  as: 


A = 


r (z0)> 

i 


r iz0) 


if  z0  * co 

if  z0  = 00 
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Where  f"  (z0)  is  the  first  derivative  of  f"  with  respect  to  z at  zQ. 

Note  that,  the  multiplier  is  same  at  all  periodic  points  of  a given  orbit.  Therefore  it  can  be  regarded  as  multiplier  of 
the  periodic  orbit. 

The  absolute  value  of  the  multiplier  is  called  the  stability  index  of  the  periodic  point.  It  is  used  to  check  the 
stability  of  periodic  points. 

Definition  2.5:  A periodic  point  Z0  is  called  attracting  periodic  point  if  | X | < 1 , and  is  repelling  if  | X | > 1 . It  is 
called  indifferent  when]  ^,1  = 1. 

Example  2.1:  Here  we  discuss  the  stability  of  periodic  points  of  the  function  f(z)  = z3  . The  fixed  points  are 

zeros  of 

/(z)  - z = 0 

Which  are  found  to  be  z}  = 0,  z2  = 1,  z3  = — 1 . Now,  f\z ) = 3z2,  so  | f\zx  )|  = 0 < 1 , therefore,  zx  = 0 
is  a stable  fixed  point.  Likewise,  z2  = 1 and  z3  = — 1 are  both  unstable  fixed  points. 

To  find  the  2-periodic  points  we  need  to  solve  the  equation; 

/2(z)  — z = 0 


The  solutions  to  this  equation  are  Zx  = 0 and  the  eighth  roots  of  unity  which  are 


1,  CO,  CO2 , CO3 , CO4 , CO5 , CO6,  CO1  where,  CO  = COS 


f 2n\  . . ( 2 

-zsm 


V o J 


V » J 


=2=(  > + ,•>. 


Since  we  have, 

(i ) f(co)=co\  f2{<o)  = f(co3)=  CO9  = CO 

(ii)  f(co2)=co6,  /2(<y2)=  /(ty6)=&>18  = co2 

(iii)  /(<y5)  = co'5  = co7,  f2 (co5 ) = f(co7)=co21  = co5 
So,  the  2-cycles  are  { CO,  CO3  j,  \a)2 , CO 6 j and  {ty5 , CO 1 } . 

As  Zj  = 0,  z2  = 1,  z3  = CO4  = — 1 are  fixed  points  of  f and  we  have  already  discussed  their  stability  we 
consider  the  points  \cok  \ k = 1,  2,  3,  5,  6,  7 }. 

Now,  / ' ( Cij k ) f ' (/  * ))  | = | (9ty8/‘  )|  = 9 > 1 . So  each  eighth  root  of  unity  is  unstable. 

After  clarifying  the  definitions  with  examples  now  we  go  back  to  the  study  of  the  maps  of  the  form 
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/C(z)  = z2+c  (1) 

For  different  values  of  the  parameter  C e C the  reason  for  choosing  the  above  map  is  that  every  quadratic 
polynomial  is  linear  conjugate  to  a map  of  this  form  and  so  it  is  representative  of  all  quadratic  polynomials.  An  important 
property  of  the  Julia  set  of  a rational  function  is  that  it  is  equal  to  the  closure  of  its  repelling  periodic  points.  A rational 

-j  where  P 

Every  polynomial  of  the  form  Cln  x"  + Cln  _\Xn  1 + ■ ■ ■ + a0  can  be  regarded  as  a rational  function.  Therefore,  the 
conceptually  easiest  way  to  define  Julia  set  is  as  follows: 

Definition  2.6:  Let  c be  any  complex  number.  The  smallest  closed  set  in  the  complex  plane  that  contains  all 
repelling  fixed  points  and  all  repelling  periodic  points  of  the  map  f ( z ) = Z2  + C is  called  Julia  set  of  the  map  f , and 
it  is  denoted  by  J c . 

For  polynomials  of  dgree  at  least  two,  oo  is  always  an  attacting  fixed  point  in  the  light  of  our  earlier  discussion  of 
stereographic  projection.  Some  authors  such  as  Goodson  [9]  use  this  fact  to  define  the  Julia  set  of  polynomials  as: 

Definition  2.7:  The  basin  of  attraction  of  co  for  the  polynomial  f (z)  having  degree  at  least  two,  is  the  open  set 

B f (go)  = jzeC: /"(z)  — > co  as  n — > ooj. 

Definition  2.8:  The  Julia  set  J(f)  of  the  polynomial  f (z)  having  dgree  at  least  2,  is  the  boundary  of  the 
open  set  B f (co)  , i.e.,  the  set  B f (go)  — B f (go). 

Definition  2.9  : The  set  K(  j ) of  all  those  points  of  C which  do  not  converge  to  go  under  iteration  of  the 
polynomial  /(z)  having  degree  at  least  two,  is  called  the  filled  in  Juia  set  of  /(z),  i.e., 

K(f)=  C-Bf(  oo). 

Note  that,  the  Julia  set  J (f)  is  also  the  boundary  d K(f)  of  the  set  K ( f ) . Thus,  for  each  z0  e j{f ) there 
is  an  open  sphere  Sr  ( Z0  j with  centre  at  z0  and  radius  r > 0 , there  is  a point  U e S r ( z0  ) such  that  iterates  of  U 

under  f converge  to  infinity  as  well  as  another  point  V e Sr  ( z0  ) such  that  iterates  of  V under  f do  not  converge  to 
infinity. 

Definition  2.10:  The  complement  of  the  Julia  set  J (f)  of  the  polynomial  /(z)  is  called  the  Fatou  set.  It  is 
denoted  by  F (/)  , i.e.,  F (/)  = C - J (/)  . 

Theorem  2.1:  For  any  polynomial,  /(z)=  anZn  + d„_lz"  1 + ■■■  + axZ  + aQ  ; an  ^ 0,  there  is  a real 


function  is  defined  as  R ( x ) = 


Q{ 


X 


( X ) and  Q ( X ) is  polynomials  with  roots  distinct  from  each  other. 
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number 

infinity. 


r such  that  if  z > r , then  f ( z ) >2  z . Furthermore,  the  iterates  of  f are  either  bounded  or  tends  to 


Proof:  Since  Cl„  ^ 0 we  choose  r = («  - i)  f T . Thus,  for  z > r we  have 

V a„  1 1 


z > ("-i) 


I I n— 1 \ ^ ^ | II  I ^ \ ^ II  I \ o I I 

z > i r =>  — \a„  z > 2 =>  — \a„  z > 2 z 


Now  we  can  make  r sufficiently  large  to  ensure  that  when  z > r , 


+ |aj|z|  + |a0|  < — | an  | | z | 


Thus 


|/(z)|  = 


II  1 n 

/ 

1 \n 

1 1 1 

1 II  1 1 \\ 

z — 

{ <*„-! 

|Z|  + an- 1 

| z | + •• 

•*  + | 1 1 Z | + 1 ££q  | 1 

^1  II  I"  ^ I II  |«  ^ 1 | ||  I”  ^ T I I 

>\an\\z\ \a„  z >—\an\\z\  >2  z 

Now,  let  us  consider  the  set  j / " (z)  : n e N | of  all  iterates  of  f . 

If  | f"  (z  )|  < r for  all  n g N then  clearly  the  iterates  of  f are  bounded. 
Otherwise,  if  | fm  (z  ) > r for  some  m g N then  we  have 

|/m+M^)|  = |/(/m+i_1(^))| 

> 2 1 fm+k~1  (z)| 

> 22 1 fm+k~2  (z  )| 


>2A'  r{z) 
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Thus,  |/'"  + A(z)|  > 2A|/"'(z)|  > 2A  r , so  |/A  (z)|— »co  as  k— >co  , i.e.  the  iterates  of  / tends  to 

infinity. 

3.  SOME  DYNAMICAL  PROPERTIES  RELATED  TO  THE  JUIA  SET 

In  this  section  we  discuss  some  dynamical  properties  of  the  Julia  sets  Jc  for  different  values  of  the  parameter  C 

Let  us  start  with  f (z)  = Z2  + C where  C is  an  arbitrary  fixed  complex  parameter. 

First  we  consider  the  case  c — 0 i.e.,  f0  (z)  = z2 . 

The  fixed  points  of f0  are  given  by 
/o  (z)  = z i.e.  z = 0,  1 

t 

As,  f0  (z)  = 2z , z = 0 is  attracting  and  z = 1 is  repelling  fixed  points  of  f0  . 

Let  z0  = re‘B , then  the  orbit  of  z0  under  f0  is  given  by: 


Thus,  the  image  of  Z0  under  fQ  is  obtained  by  squaring  the  magnitude  of  Z0  and  doubling  the  angle  (argument  ) 
of  Z0  . Therefore,  the  behavior  of  the  orbit  of  Z0  depends  on  its  magnitude. 

If  r < 1 , then  r2  — > 0 as  n — > go  , therefore,  the  orbit  of  z0  will  tend  to  the  origin  which  is  an  attracting  fixed 
point  of f0. 

Ifr  = 1 , that  is,  if  z0  is  on  the  unit  circle  | z | = 1 , then r 2 = 1 , so  the  orbit  of  z0  will  remain  on  the  unit  circle 
and  on  each  iteration  f0  will  doubles  the  angle  (argument)  of  z0 . 

Ifr  > 1 , then  r — >oo  as n — » go , and  so  the  orbit  of  z0  will  tends  toco,  so  the  orbit  in  this  case  will  be 

unbounded.  This  situation  has  been  shown  in  the  following  cob-web  diagram. 


www.iaset.us 


editor@iaset.us 


108 


Arun  Mahanta,  Hemanta  Kr.  Sarmah,  Ranu  Paul  & Gautam  Choudhury 


Figure  3.1:  The  Cob-Web  Diagram  for  /0 

Further,  the  following  Table  1 gives  the  iterated  values  of  z0  under  f0  in  each  of  the  above  mentioned  situations. 


Table  3.1 


z0 

Z,  = f0(z„) 

Z2  = f0(zi) 

Z3  = f0(z2) 

Z4=f„(z3) 

Z5=f0(z4) 

| z | < 1 

Modulus 
= 0.8 

Argument 
= 10° 

Modulus 
= 0.64 
Argument 
= 20° 

Modulus 
= 0.4096 
Argument 
= 40° 

Modulus 
= 0.1678 
Argument 
= 80° 

Modulus 
= 0.0281 
Argument 
= 160° 

Modulus 
= 0.0008 
Argument 
= 320° 

I z | = 1 

Modulus 

= 1 

Argument 
= 10° 

Modulus 

= 1 

Argument 
= 20° 

Modulus 

= 1 

Argument 
= 40° 

Modulus 

= 1 

Argument 
= 80° 

Modulus 

= 1 

Argument 
= 160° 

Modulus 

= 1 

Argument 
= 320° 

I z | > 1 

Modulus 
r = 1.5 
Argument 
0 = 50° 

Modulus 
= 2.25 
Argument 
= 100° 

Modulus 
= 5.06 
Argument 
= 200° 

Modulus 
= 25.6 
Argument 
= 40° 

Modulus 
= 655.36 
Argument 
= 80° 

Modulus 
=431440 
Argument 
= 160° 

The  iteration  of  three  initial  points  mentioned  in  the  Table-1  is  graphically  represented  in  figure -4. 


Figure  3.2:  Iteration  of  the  Point  z0  under  f0  as  Shown  in  Table-3.1 
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Since  the  doubling  map  on  the  unit  circle  is  chaotic  [4],  f0  is  chaotic  on  the  unit  circle.  This  circle  is  the 


itineraries  that  go  to  infinity.  This  separating  curve  (circle)  is  the  Julia  set  of  f0  and  clearly  the  interesting  dynamics  of 
f0  takes  place  on  it.  At  this  point  it  is  noteworthy  that  although  fQ  exhibits  chaotic  behaviour  on  its  Julia  set,  its  study 
may  lead  one  to  believe  that  Julia  sets  are  nice  smooth  curves.  In  fact,  a variation  in  the  parameter  C give  rise  to  quadratic 
maps  whose  Julia  sets  are  fractals.  Ruelle  [17]  has  shown  that  for  fc  with  small  |c  |,  the  fractal  dimension  of  its 
corresponding  Julia  set  is  approximately 


And  so  is  indeed  a fractal  in  these  cases. 

Next,  it  is  to  be  noted  that  the  fixed  points  of  fc  are  given  by 

1 ± V 1-4 c 

z = 

2 

11 

Clearly,  number  of  fixed  points  will  differ  depending  on  c — — or  c A — 

4 4 

Now,  we  discuss  the  nature  of  attracting  and  repelling  fixed  points  of  ( 1 ) through  the  following  theorems: 
Theorem  3.1:  fx  Have  neither  attracting  nor  repelling  fixed  points 


boundary  between  the  set  of  initial  points  in  the  complex  plane  with  itineraries  that  approaches  the  origin  and  those  with 


d n — 1 + 


I I2 

I c I 


C 


4 log  2 


4 


Proof:  Fixed  points  of  fx  is  given  by  the  quadratic  equation: 

4 


/l  0)  - z 


4 


Thus,  the  only  fixed  point  of  fx  is  at  z = — . 

4 2 


N°w>  f\  - 

4 v 2; 


= 1 , which  shows  that  z = — is  indifferent  fixed  point  of  f t . 

2 4 


Therefore,  fx  has  neither  attracting  nor  repelling  fixed  point.  ■ 
4 


Theorem  3.2:  Ifc  ^ — , then  at  least  one  fixed  point  of  f is  repelling. 
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Proof:  Fixed  points  of  f is  given  by  the  quadratic  equation: 
fc{z)  = z =>z2  + c = z 

Thus,  the  fixed  points  of  (1)  are  the  roots  Zx  and  Z2  of  the  quadratic  equation 

z2  — z + c = 0 And  therefore,  z,  — — + — ,/l  — 4c  andz,  = — -A/l  — 4 c . 

22^  22^ 

1 , 

From  this  it  is  clear  that  for  c ^ — , J c has  two  distinct  fixed  points  viz.,  zl  and  Z2  . 

r 

Now  fc  (z)  = 2 z , therefore, 

2,  = fc  (zj)  = 1 + i/l-4c  = 1 + //,  and2,2  = /c  (z2)  = 1 - ^1-4 c = 1 - JU, 

Where  /U  = -yjl  — 4c 

As  c ^ — , u 0 . 

4 

Consider  the  unit  circle 
| z | = 2 •••  (2) 

Two  cases  might  arise  in  case  of  position  of  a point  with  respect  to  the  circle  (2)  depending  on  the  value  of  /U 
which  are  discussed  below: 

Case  1:  When  /u  lies  outside  the  circle  (2). 

In  this  case  both  the  distances  of  /U  and  — jU  from  1 are  greater  than  unity,  i.e. 

1 1 - /U  | > 1 And  1 1 + /U  | > 1 
=>  | X2  | > 1 And  | Al  | >1 

Therefore,  both  Z,  and  Z2  are  repelling  fixed  point  of  (1) 

Case  2:  When  ju  lies  inside  or  on  the  circle  (2). 

In  this  case,  as  /d  ^ 0 at  least  one  of  the  distances  of  /d  or  — fd  from  1 is  greater  than  the  unity,  i.e.  at  least  one 
of  | 1 — fd  | or  | 1 + fd  | is  greater  than  one.  Therefore,  at  least  one  of  Z,  or  Z2  is  repelling  fixed  point  of  (1).  ■ 

Theorem  3.3:  f\  (z)  = z2  4 has  repelling  2-cycles. 
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Proof:  It  is  already  shown  in  theorem  3.1  that  the  only  fixed  point  of  f]  (z)  is  z = — which  is  neither  attracting 


nor  repelling. 


Now  the  2-cycles  of  fl  are  given  by 

4 


fl  0)  = z 

4 


f „ O 
z2  + — 

V 4y 


— = z =>  z4  + —z2  — z + — = 0 (3) 


4 


16 


1 2 

Since  — is  fixed  point  of  fx  (z)  , it  is  also  a fixed  point  of  fl  " (z)  and  therefore, 

^ 4 4 


f 1 > 

Z : 

v 2y 


is  a factor  of  left 


hand  side  of  (3).  On  factorization  we  get  from  (3); 

( iY 

z 

V 2 J\  Z H OJ 


3 1 2 3 5^ 

z + — z + — z 

2 4 8 


= 0 


f 


It  turns  out  that 


1 


z — 


■J 


3 1 2 3 5 1 

are  again  a factor  of  z -I z H z , i.e.  — is  a double  root  of  the 

2 4 8 2 


equation  (3),  which  leads  to  the  following  factorization: 

- 0 


f 0 

2 

f 2 5) 

l 2; 

l 4 ) 

11  . 2,  S 1 1 1 

By  solving,  z = — , ± I Therefore,  the  fixed  points  of  /,  (z)  are  — , 1-  I,  and I . But, 

2*  2*  , 2*  2t  2* 


1 1 P / \ 1.1., 

neither 1 hi  nor i is  a fixed  point  of  /,  (z) , this  leads  us  to  conclude  that  < hi, l>  is  a 2- 

2 2 ^ l 2 2 

cycle  of/,  (z). 


To  show  this  2-cycle  is  repelling  we  first  find  the  multiplier  X and  then  use  the  definition  1.5  as  follows: 


\x\  = 

f \ 

A2 

( 1 

b l 

— 

r 

A 

f 

A 

r i 

b l 

^ , 
f\ 

f 1 

b l 

V 4 ) 

V 2 ) 

4 

V 4 

\ 2 j 

/ 4 

K 2 ) 

A 

4 V 


I 

2 


A 

J 4 


b l 


2 

1 

1 

.2 

r i 

— + 1 

V 2 J 

V 2 ) 

; Since,  / 

i 

; As,  /,  (z)  = 2z 


( 1 

b l 

l 2 J 


1 

i 

2 
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= |(-l-2/)(-l  + 2/)|  = 5 

As  | A,  | > 1,  the  2-cycle  | — ZJ  is  repelling.  ■ 

4.  SOME  PROPERTIES  RELATED  TO  THE  STRUCTURE  OF  JULIA  SETS 

In  this  section  we  will  prove  some  results  relating  to  the  structure  of  Julia  sets. 
Property  4.1:  For  every  complex  number  C , the  Julia  set  of  f is  non-empty. 


Proof:  By  theorem  3.2,  fore  ^ , Jc  has  a repelling  fixed  point  and  by  theorem  3.3,  fore  = — ,JC  has  a 

repelling  period-2  cycle.  As,  J (fc)  i.e.  Julia  set  of  fc  is  the  closer  of  the  set  of  all  repelling  periodic  points  of  fc 
hence  j(fc  W-- 

Property  4.2:  If  | z | > | c | + 1 , then  the  orbit  of  z for  fc  is  unbounded. 

Proof:  When  c = 0 , the  iterates  of  z are  just  positive  powers  of  z as  f0  ( z ) = z2 , which  is  clearly  unbounded 
for  I z I > 1 . 


Suppose,  c 0 and  z > c + 1 , then 


\fc{Z)  I = \z2  + C = \z\ 

c 

Z H 

> I z | 

1 1 c 
1 z 1 

z 

z 

> 


z 


c 

c 

+ 1 

Where,  r — 


c | 

cUT 


> 1 as  c 


1 

c 

1 

c 

+ 1J 

1 

c 

c 

+ 1 

It  follows  that  fcn(z)  > r"  | z | — > oo  as  n increases  without  bound.  Therefore  if  | z | > | c | + 1 , then  the 
iterates  of  z form  an  increasing,  unbounded  sequence,  so  the  orbit  of  z is  unbounded.  ■ 

Definition  4.1:  A set  S is  called  invariant  under  the  function  f if  i.e.  the  image  of  S under  f 

is  contained  in  S . 

Property  4.3:  The  Julia  set  J of  any  complex  polynomial  / is  completely  invariant  under/.  That  is, 
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Proof:  Suppose, 

zeJ= i>  fk+l(z)  = fk  (/  (z)  ) Not  tends  to  infinity  as  k — > go 

Recall  that,/  is  the  boundary  of  the  basin  of  attraction  of  infinity  and  z is  a repelling  point  of  f . 

As  f is  continuous  at  z , we  can  find  Un  — > z such  that 

fk+l  (u„  ) = fk  (f  ( Un ))  Tends  to  infinity  as  k — > co  for  all/; 

Therefore,  f (z)  is  not  on  the  basin  of  attraction  of  infinity  and  we  can  find  points  f ( Un  )as  close  to  f (z) 
as  we  wish  such  that  f (un  ) are  on  the  basin  of  attraction  of  infinity.  From  which  one  can  conclude  that 

/(z)gJ  (4.3.1) 

This  implies  that 

(4.3.2) 

Next  consider  the  pre -image  z0  of  z i.e.  f (z0  ) = Z . By  mapping  properties  of  / , we  may  find  v„  — > z0  with 
f(vn)  = un. 

fk  (z0  ) = /*“'  (z)  Not  tends  to  infinity  as  A:  — » oo 
And  fk  (vn)  = fk  1(un)  tends  to  infinity  as  k — » oo  for  all/; 


Hence  z0  eJ^>z  = / ( z0  ) e / [ / ] (4.3.3) 

J ^ /[J]  (4.3.4) 

(4.3.2)  and  (4.3.4)  implies  that 

J = /[J]  (4.3.5) 

From  (4.3.1) 

z G /"*[*/]  And  therefore,  J c (4.3.6) 

From  (4.3.3) 

/_1(z)  e J And  therefore,  c J (4.3.7) 

(4.3.6)  and  (4.3.7)  implies  that 

J = /-‘[J]  (4.3.8) 


(4.3.5)  and  (4.3.8)  together  constitute  the  result.  ■ 
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Property  4.4:  The  orbit  of  z e J ( fc  ) is  bounded. 

Proof:  Letz  e J ( /'.  ).  Now  if  z is  a periodic  point  of  f , then  necessarily  the  orbit  of  z is  bounded  and 
therefore  by  property  4.2,  | Z | < | C | + 1 . Again  J ( /'.  ) is  the  smallest  closed  set  containing  all  repelling  periodic  points 
of  fc  , therefore,  any  z e J \fc  ) also  has  the  property  that  | Z | < | C | + 1 . Now  by  property  4.3,  the  iterates  of  j(fc  ) 
are  precisely  J (fc  ) • Hence  the  iterates  of  z are  also  bounded  by  | c | + 1 . ■ 

Proposition  4.5:  The  Julia  set  J c of  the  function  fc  is  compact  for  allc  e C . 

Proof:  By  theorem  1.1,  there  exist  a real  number  r such  that  for|  z | > r , 

|/c(z)|  ^ 2z 

Thus,  for  | z | > r , f"  (z)  — » oo  as  n — » oo  Since,  Julia  set  J c of  fc  is  the  boundary  of  the  set  of  those 
points  of  C which  do  not  converge  to  co  under  iteration  of  f , so  we  must  have 

J c C Sr  ( 0 ) I.e.  the  sphere  of  radius  r and  centre  at  origin 

=>  | z | < r For  all  z e J c 

Which  shows  that  J c is  bounded. 

Next  we  need  to  show  that  J c is  closed. 

Let  z e5  ( °o  ),  the  basin  of  attraction  of  co  under  fc . Therefore,  there  exist  k e N such  that  fck  (z  ) > r . 
Thus,  to  each  z G B ( oo  ) , there  exist  an  open  sphere  S s ( z ) such  that 

ze5£(z)  =>  fk  ( z ) — » oo  As  A:— >oo 

=>  z e B (oo ) 

Therefore,  z G S£  (z  )c  B [ oo  ) which  shows  that  B ( oo  ) is  an  open  set. 

Now  asi?(co)  = C — K ( fc  ) , so  the  filled  in  Julia  set  K(fc  ) is  closed.  Again  the  boundary  of  K ( fc  ) 
is  the  Julia  set  J c we  have  that  J c is  closed. 

Hence  J c is  compact.  ■ 

Proposition  4.6:  J c is  symmetric  about  the  origin  for  all  C G C . 
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Proof:  Suppose  zo  6 Jt.  By  invariant  property  of  Julia  set  it  is  clear  that 

/c(Zo)  e Jc  =>  Z0  + C G J c ^ f c (Z0  ^C)  £ ^ {Z0>  _Zo}£  ^ _Z0  6 >^c 

Hence  the  Julia  set  J c of  fc  is  symmetric  about  the  origin.  ■ 

Proposition  4.7:  For  any  real  values  of  c , the  Julia  set  J is  symmetric  about  both  real  and  imaginary  axes. 
Proof:  Recall  that  any  set  iS  (Z  C is  symmetric  about  the  real  axis  if  the  conjugate  of  each  element  of  S 

belongs  to  the  set  S i.e.  z e S =>  z e S . We  first  show  that  f"  ( z ) = f”  (z ) . 

Let  Zj  = flj  + ib j 

Now,  fc  (zj  )-a2  - b j2  +c  +2albli  — a2+ib2  (say) 
fc  (z)  = «12  - V +c  -2a!  bx  i = a2  -ib2  =fc  (z1 ) 


Suppose  fc'~'(z1  ) = /t.'“1(z1  ) and  //“‘(Zj)  = ar  +ibr 


f \ 

f \ 

Now,  fj  ( z,  ) = fc 

*i) 

= f. 

//■*  (zl) 

V V 

L J 

= a,2  -b2  +c-2ar  br  i = //  (z, ) 


Thus,  by  induction, 

f \ 


i:  z = fcn(z) 


For  all  12  G iV 


V J 

Now  suppose, 

z e Kc , the  filled  in  Julia  set  of  f 

=>  k"(-)L°°  Isbounded 

=>  j /c”  (z)  r °°  Is  bounded 


ff(z)\  °°  Isbounded  =>  z e Kc 


Thus  Kc  and  hence  its  boundary  J c is  also  symmetric  about  the  real  axis. 

Symmetry  of  J c about  the  imaginary  axis  immediately  follows  from  the  fact  that  it  is  symmetric  about  the  real 
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axis  as  well  as  the  origin.  ■ 

5.  VISUALIZING  JULIA  SETS 

To  see  how  one  can  plot  Julia  sets,  one  will  look  at  the  simplest  method  called  "The  escape  criterion  method", 
which  generates  the  filled  in  Julia  sets.  Julia  sets  do  not  have  interiors;  however,  this  method  often  ends  up  plotting  the 
interior  along  with  the  Julia  set.  This  method  is  based  on  the  following  facts: 

At  each  step  of  iteration  of  the  map  f (z)  = z~  + C , geometrically  we  square  the  modulus  and  double  the  angle 
of  the  number  to  which  it  operates  and  then  shift  by  the  parameter  C . If  C is  small  the  squaring  part  dominates  the  shift  so 
the  behavior  is  not  very  different  from  the  case  c — 0 , discussed  above.  Just  as  this  case,  two  types  of  behavior  are 
possible  depending  on  the  starting  point  (seed).  The  first  type  of  points  lie  in  the  basin  of  attraction  of  the  fixed  points  or 
the  attracting  cycles  and  the  orbits  of  the  second  type  tend  to  infinity.  The  Julia  set  separates  these  two  very  different  types 
of  behavior.  If  a point  is  in  the  Julia  set,  there  are  arbitrarily  close  points  that  iterates  to  infinity  but  also  arbitrarily  close 
points  that  do  not  wander  far  off  under  iteration.  So,  our  first  method  for  plotting  the  Julia  set  is  to  color  points  with  a color 
if  their  orbits  tend  to  infinity  and  to  color  them  by  another  one  ( contrasting  to  the  first  color)  if  they  do  not.  The  boundary 
between  these  two  regions  is  then  the  Julia  set. 

Now,  our  problem  is  converted  to  find  out  those  points  whose  orbits  tend  to  infinity  which  is  tackled  by  the 
following  theorem: 

Theorem  5.1:  If  | C | < 2,  Then  the  orbit  of  the  points  lie  outside  the  circle  of  radius  2 i.e.  the  set  of  points 
| z:  | z | > 2},  escape  to  infinity. 

Proof:  Let  ii  = {z  : | z | > 2j  . Then  for  any  z e E we  have, 

| z | > | c | 

Now,  | z2  | = | c + z2  — c|  < | c | + | z2  — c | 

=>|z2|  — c | < | z2  — c | =>|z2|  — | z | < | z2  — c | 

Replacing  C by  — C, 

(I  z | — i)|  z | C | z2  + c | = | fc  (z)  | 

Since  | z | > 2 , therefore,  | z | — 1 > 1 and  hence  one  can  write  | z | — 1 = 1 + 8 where  5 > 0 . 

M < (1  + s)\z\  < I fc (X>  I 

Thus,  z will  move  further  from  origin  by  the  action  of  f 


Impact  Factor  (JCC):  2.0346 


NAAS  Rating  3.19 


Julia  Set  and  Some  of  its  Properties 


117 


Again, 


fAfAA) I > (i  + <s)|A(z)|  > (i  + <5)2U 


Therefore,  f (z)  is  even  further  from  the  origin  than  f (z)  . 
Continuing  in  this  fashion,  one  can  find 

fc  (Z1  > (!  + <?)"  I z I 


Recall  that  5 > 0 , so  thatl  + 8 > 1 . Hence  the  number  (l  + 8 )"  — »oo  as  n — > co  so  it  follows  that 


//  M 


— > oo  as  n — >oo  i.e.  the  orbit  of  z tend  to  infinity,  i 


In  view  of  the  above  theorem  the  algorithms  that  produce  the  Julia  sets  shown  in  the  figure  5. 1 and  5.2  is  as 

follows: 

Step  1: 

Input  Cj  and  C2  where  C = Cl  + ic2 

Step  2: 

Select  a 200  x 200  grid  in  the  plane. 

Step  3: 

For  each  point  Z0  in  this  grid,  compute  the  first  20  points  on  the  orbit  ofz0 . Check  at  each  point  of  the  iteration 
whether  the  corresponding  iterated  value  lies  outside  the  circle  of  radius  2 . 

Step  4: 

If  any  point  on  the  orbit  lies  outside  the  circle  of  radius  2 , then  stop  iteration  and  color  the  point  Z0  by  purple 

color. 

Step  5: 

If  all  20  points  on  the  orbit  lie  outside  the  circle  of  radius  2 , then  keep  the  original  point  z0  without  any  color 
(i.e.  keep  it  white) 
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(a) 


(b) 


(c) 


-1  0 1 2-2-1  0 1 2-2-1  0 1 


(d) 


(e) 


(f) 
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(a)  - 0.1,  (b)  0.25,  (c)0.4,  {d)  0.65,  (e)0.8,  (/)  0.9,  (g)-0.5,  ( h)-0J , (0-1) 


G) 


(k) 


(1) 


(m) 


(n) 


(o) 


(p)  (q)  (0 

Figure  5.2:  Filled  in  Julia  Set  for  the  Parameter  (c)  Values 

O') -1.25,  (k)  -1.5,  (0  -1.9,  (rn)  -0.3  i,  (;?)  -0.5/,  (o)  /,  ( p ) 0.1  + 0.5/,(?)  -0.1  + 0.75/,  (r)- 0.1 0.9/). 
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Though  this  method  gives  beautiful  pictures  of  Julia  sets,  it  actually  gives  the  interior  together  with  the  Julia  sets 
i.e.  filled  in  Julia  sets  instead  of  showing  the  actual  Julia  sets  alone.  Another  disadvantage  of  the  algorithm  is  that  it  takes 
more  computational  time.  To  overcome  these  disadvantages,  there  are  various  computation  methods  for  visualizing  Julia 
sets  in  computer  screen  such  as  Boundary  Scanning  Method  (BSM),  Inverse  Iteration  Method  (IIM),  and  Modified  Inverse 
Iteration  Method  (MUM)  etc.  These  methods  are  explained  and  compared  in  [8],  [19]  and  [3].  Among  these  MUM  is  the 
best  method  for  both  computation  time  and  remarkable  picture  quality.  This  method  is  discussed  in  detailed  by  Peitgen  and 
Richter  [15].  The  algorithm  for  this  method  is  as  follows: 

Step  1: 

Choose  any  real  number  c such  that  — 2 < c < 2 . 

Step  2: 

Find  any  repelling  point  x0  of  fc  as  the  starting  point  and  set  S0  = { x0  } . 

Step  3: 

Choose  the  number  of  iteration,  say  n . 

Step  4: 

Find  the  image  of  S0  under  and  g2  where gj  (z)  = Jz  — c , and  g2  (z)  = --Jz  - c 

Step  5: 

SetS,  =g,[S0]ug2[S0]. 

Step  6: 

Find  the  image  of  51!  - S0  under  gj  and  g2  and  setS2  = -50]  u g,[5(  -S0]. 

Step  7: 

Repeat  the  process  for  n times  and  finally  get  Sn=  g(  [ 5K_1  - Sn_2  ] u g2  [ Sn_l  - Sn_2  ] 

Step  8: 

fl 

The  set  St  will  give  the  Julia  set  for  fc  approximately. 

;=o 

Using  this  algorithm  Julia  sets  for  the  parameter  value  c=— 1,  c = 0. 123  + 0. 745/,  c=—  0.5+  0.55/, 
c = 0.4+  0.3/  And  c = i are  drawn  which  are  shown  in  the  figure  5.3  below: 
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(e) 

Figure  5.3:  Julia  Sets  for  the  Parameter  (c)  Values 

(a)  — 1 , (b)  0.123  + 0.7451,  (c)  -0.5  + 0.551,  (d)  0.4  + 0.31  (e)  1 


6.  CONCLUSIONS 

For  all  parameter  values  c such  that  | C | < 0.25 , effect  of  shift  by  c does  not  change  the  picture  of  the 

respective  Julia  set  much.  In  this  case  the  Julia  set  is  a quasi  circle  or  a loop.  As  soon  as  the  magnitude  of  c increases  the 
shifting  part  of  the  function  become  significant,  its  picture  becomes  strange.  For  example,  the  Julia  set  for  the  parameter 
value  C = — 1 consists  of  a central  loop  surrounded  by  many  smaller  touching  loops  which  in  turn  are  surrounded  by  even 
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smaller  loops,  and  so  on,  see  figure  5.3(a).  Moreover,  the  Julia  set  for  the  parameter  value  c close  to  — 1 has  a similar 
structure  of  touching  loops  i.e.  the  configuration  is  stable  for  C = — 1 . The  sets  for  other  parameter  values  yields  more 
surprises.  For  c close  to  —0.123  +0.745/  the  Julia  set  is  again  formed  by  a hierarchy  of  loops,  but  now  they  meet  three 
at  a point  as  shown  in  figure  5.3(b).  For  this  case  fc  has  an  attracting  period-3  orbit.  Further,  increasing  the  value  of  c 
nearer  to -0.5  + 0.55/,  the  loops  meet  five  at  a point,  figure  5.3(c).  Again,  choosing  the  value  of  c as  — /,  the  loops 
collapse  and  the  Julia  set  takes  a dendrite  or  'twig-like'  shape,  see  figure  5.3(e).  As  the  value  of  c becomes  larger  a 
dramatic  change  in  the  nature  of  the  Julia  set  takes  place.  For  example,  when  c = 0.4  + 0.3/  the  connectedness  of  the 
Julia  set  completely  breaks  up  and  it  become  a totally  disconnected  one,  figure  8(d).  In  this  situation,  the  Julia  set  becomes 
a fractal  dust  and  is  identical  to  the  filled  in  Julia  set.  Any  point  not  in  the  Julia  set  of  this  type  iterates  to  infinity,  as  no 
other  destination  is  possible. 
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